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A  Perturbation  Theory  for  Light  Diffraction  from 
a  Bigrating  with  Multiple  Surf ace-Polariton  Excitation 


N.  E.  Glass 
Department  of  Physics 
Naval  Postgraduate  School 
Monterey,  CA  93943 


Abstract 

A  perturbation  theory  has  recently  been  developed  for 
treating  the  diffration  of  light,  with  surface-polariton  resonant 
excitation,  from  a  bigrating  surface.  The  theory  is  an 
approximation,  to  first  order  in  the  grating  height,  of  an  earlier 
theory  based  on  the  Rayleigh  hypothesis.  The  perturbation  theory 
holds  for  arbitrary  polarization  and  for  arbitrary  plane  of 
incidence  with  respect  to  the  grating  directions.  It  is,  however, 
limited  to  treating  only  the  possibility  of  two  simultaneous 
resonant  evanescent  waves,  whereas  four  are  possible  at  normal 
incidence  on  a  square  grating,  and  three  are  possible  on  a 
rectangular  grating.  The  present  work  generalizes  the  earlier 
work  to  allow  for  a  four-fold  (or  three-fold)  resonance.  This 
extension  also  allows  one  now  to  determine  the  complex  dispersion 
relation  for  surface-polaritons,  for  wavevectors  at  the 
intersection  of  two  Brillouin  zone  boundary  lines  (not  possible  in 
the  previous  formulation) . 
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I.  Introduction 


The  study  of  electromagnetic  wave  scattering  from  corrugated 

metal  surfaces,  with  resonant  excitation  of  surface-plasmon- 

polaritons,  has  many  important  technological  applications.  One  of 

these  is  the  creation  of  electromagnetic  radiation  from  an 

electron  beam  passing  above  a  grating  surface  --  the  Smith-Purcell 

effect'  .  In  many  applications,  the  enormous  enhancements  in  the 

absorptance  at  the  surface  and  in  the  field  intensity  in  the 

surface  region,  due  to  the  grating  coupling  of  an  incident  wave 

into  a  surface-polariton  (or  surface-electromagnetic-wave,  SEW) , 

can  be  further  enhanced  by  use  of  a  bigrating  (a  surface  periodic 

in  two  different  directions)  in  place  of  a  classical  grating 

(periodic  in  one  direction) .  It  was  pointed  out  in  some  of  the 

f2-3) 

earlier  theoretical  studies  on  bigratings'  '  that  the  resonant 

absorption  of  unpolarized  light  (e.g.,  of  sun-light  in  solar 

absorbers)  will  be  greater  on  a  bigrating  than  on  a  classical 

grating.  When  the  plane  of  incidence  is  normal  to  the  groove 

direction,  only  the  p-polarized  component  of  the  incident  wave  can 

excite  the  surface-polariton  (which  is  itself  a  p-polarized 

surface-wave) ;  the  s-component  does  not  couple  to  the  surface-wave 

at  all.  By  contrast,  the  bigrating  can  couple  both  the  p-  and  the 

s-  components  of  the  incident  wave,  regardless  of  the  orientation 

of  the  plane  of  incidence.  Moreover,  as  demonstrated 

f  4  ^ 

experimentally  by  Inagaki  et  al.  ,  the  bigrating  can  couple  the 
incident  wave  simultaneously  into  two  surface-polaritons ,  leading 
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to  an  extra  absorptance,  for  incidence  orientations  in  which  the 
classical  grating  would  only  excite  a  single  polariton.  These 
same  processes  indicate  mechanisms  by  which  the  Smith-Purcell 
radiation  can  also  be  enhanced  by  using  a  bigrating  in  place  of  a 
classical  grating. 

The  first  theory  of  the  Smith-Purcell  effect  that  was  both 

mathematically  exact  and  physically  realistic  was  the  work  of  P.M. 
CS)  . 

van  den  Berg'  .  This  theory  used  a  Green's  function  approach, 

( S) 

proposed  by  di  Francia^  ,  in  which  the  component  evanescent 

waves,  which  make  up  the  Fourier  expansion  of  the  field  of  the 

electron  current,  are  seen  as  exciting  the  various  diffracted 

orders  of  the  classical  grating. 

One  shortcoming  of  the  van  den  Berg  theory  is  that  it 

considers  the  metal  on  which  the  grating  is  ruled  to  be  a  perfect 

conductor.  This  approximation  precludes  the  possibility  of 

accounting  for  the  effects  of  the  surface-polaritons  and  their 

( 1 ) 

resonant  excitation.  Recently,  Chuang  and  Kong'  '  developed  a 
theory  which  considers  a  sinusoidal  grating  on  silver,  described 
by  a  complex  dielectric  function,  and  uses  an  extended  boundary 
condition  method  to  calculate  the  radiated  energy. 

Similar  techniques  are  in  principle  applicable  to  the 
treatment  of  Smith-Purcell  radiation  from  a  bigrating,  but  the 
numerical  difficulties  can  become  prohibitive:  the  computing  time 
necessary  to  calculate  the  diffracted  field,  arising  from  each 
Fourier  component  in  the  incident  field,  will  be  much  greater  for 
the  two-dimensional  case,  the  bigrating,  than  for  the 
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one-dimensional  case,  the  classical  grating.  Thus  perturbative 
techniques  are  desired. 

A  non-perturbative  theory  of  light  diffraction  from  a 

bigrating,  with  sur face-polariton  resonances,  was  developed  by 

f  8 ) 

Glass,  Maradudin,  and  Celli'  .  This  theory  uses  the  vectorial 
form  of  Kirchhoff's  integral  to  eliminate  the  field  in  the  metal 
and  thereby  to  cut  in  half  the  size  of  the  matrix  equation  to  be 
solved.  There  are  nevertheless  very  large  matrix  equations  to  be 

solved  with  this  method:  from  98x98  up  to  243x243  in  Ref.  8. 

.  ( 9 1 

Thus,  a  perturbation  theory  was  recently  developed  by  Glass'^  , 

which  applies  the  coupled-mode  theory  of  Glass,  Weber,  and 

Mills^^^^  to  the  exact  theory  of  Glass,  Maradudin,  and  Celli^®^. 

This  perturbation  theory  reduces  the  computations  to  a  4x4  matrix 

equation,  by  retaining,  to  first-order,  the  two  amplitudes  (p-and 

s-  polarized  components)  of  the  specularly  diffracted  beam  and  two 

evanescent-wave  amplitudes  (p-components  only)  corresponding  to 

the  possible  simultaneous  resonant  excitation  of  two  surface- 

polaritons.  These  two  surface-polaritons  may  be  propagating  in 

non-collinear  directions  (unlike  the  case  for  the  classical 

grating,  in  which  two  polar itons  can  be  excited  only  for  normal 

incidence,  and  then  are  propagating  in  exactly  opposing 

directions).  This  theory  as  implemented  in  Ref. 9  has  one 

shortcoming:  it  cannot  treat  the  case  of  more  than  two 

surface-polaritons  excited  simultaneously.  For  example,  it  cannot 

treat  the  cases  of  normal  incidence  on  a  square-grating  when  there 

are  four  surface-polaritons  excited  simultaneously  or  of  oblique 
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incidence  on  a  rectangular  grating  when  there  are  three 
simultaneous  excitations.  This  shortcoming  could  be  an  important 
one  if  one  wishes  to  use  these  perturbative  approaches  in 
calculating  the  Smith-Purcell  radiation.  Fortunately,  it  is  not  a 
difficult  problem  to  correct.  A  straight-forward  extension  of  the 
theory  of  Ref.  9  will  produce  a  6x6  matrix  equation  for  four 
resonant -wave  (p-)  amplitudes  and  the  two  specular  beam 
amplitudes.  It  is  the  purpose  of  this  report  to  present  the 
results  of  the  extension  of  the  perturbation  theory  of  Ref.  9  to 
treat  the  cases  in  which  there  are  up  to  four  surface  excitations. 

II.  Theory 

The  physical  geometry  is  as  follows.  The  nominal  flat 

surface  is  in  the  plane  X2=0.  The  surface  profile  is  described 

by  X2=C(X||)  [where  x^^=x^x^+x^x^  ].  The  region  X2>f(X||)  is 

vacuum;  X2<C(X||)  is  dielectric,  described  by  the  complex 

dielectric  function  e  (w)  =ej^(w) +ie ^  (w)  .  The  surface  profile 

function  is  periodic  in  two  directons.  Light  of  frequency  w  and 

wavevector  Ic  is  incident  from  the  vacuum  (with  angle  of 
•  •  •  • 

incidence  0).  The  preelection  of  k  onto  the  surface  is  k|| , 

A 

which  makes  an  angle  «#>  with  x^. 

The  electric  field  in  the  vacuum  is 
f(x)  =  exp[iiC|| ‘Xu  -  ia^(w,k||  )X2]  + 

+  2  exp[iSg-X||  +  ia^((j,K^)X2]  .  (1) 

In  the  above  equation  we  have: 


5 


(2) 


%  =  Sii  +  3  , 

where  G  is  a  reciprocal  lattice  vector  defined  by  the  bigrating 


periodicity,  which  for  a  simple  square-lattice  grating  is 


2tt  ^  ^ 


for  mj=0,±l,±2 , . . . 


("a"  is  the  period);  and  where 

(jl)  ^  ^ 

k||  =  Kq  =  -  {cos(p  +  sin<t>  X2)sin0 
Also  appearing  in  Eq.(l)  is 


ao(a.,K^) 


,2,2  x,,2,l/2 

(o  /c  -  K^)  ^  , 

.  ,^2  2  ,  2,1/2 

l(Kg  -  w  /c  )  '  , 


for  <  (o^/c^ 

Cj 

2  2  2 
for  >  w^/c 

Cj 


(3) 

(4) 

(5a) 

(5b) 


The  amplitudes  of  the  incident  and  scattered  fields  are: 

3  a  (!"k„)  *  (^3  >=11)®! 


^^(w,)c  )  =  [k„+x 


and 


(6) 


r(o,%)  =  [K^  -  X3  ]A„(o,%)  +  (X3  X  Kg)A_L(c.,%)  (7) 

The  p-  and  s-  components  for  the  incident  wave  are  expressed  in 
terms  of  B||  and  ,  and  for  the  scattered  waves  in  terms  of 

the  Rayleigh  coefficients  A||  and  Aj^ 

The  Rayleigh  coefficients  are  the  unknowns  in  the  problem. 

To  find  them  one  could  expand  the  field  in  the  metal  in  a  Rayleigh 
expansion  and  then  apply  the  boundary  conditions  at  the  grating 
surface  (this  is  the  Rayleigh  method) .  The  method  used  in  Ref.  8 
reduces  the  number  of  equations  in  half  by  employing  the  vectorial 
equivalent  of  the  Kirchhoff  integral  to  eliminate  the  field  in  the 
metal.  The  effective  boundary  condition  that  results  then  yields 
two  infinite  sets  of  linear  algebraic  equations  for  the  unknown 
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coefficients.  These  are  Eqs.  (2.18a)  and  (2.18b)  in  Ref.  8  and 
equations  (8a)  and  (8b)  in  Ref  9. 

The  equations  for  the  Rayleigh  coefficients  contain  the 
integral : 

I(“|G)  =  -^  J'J'  exp[-iG-X||-iaC  (x,| )  ]  ,  (8) 

c  cell 

where  a^  is  the  area  of  a  unit  cell  on  the  bigrating  surface,  over 
which  the  integral  is  performed. 

The  perturbation  theory  begins  by  expanding  the  integrand  of 
"I"  to  first  order  in  C,  thus  yielding: 

I(q:|G)  =  -  iaHG)  (9) 

where 


C(^)  =  d^X||C  (X||  )exp(-iG-X|| )  .  (10) 
c  cell 

Under  the  assumption  that  the  dimensionless  parameter  af  is  small, 
we  can  substitute  equation  (9)  into  the  equations  for  the  A's,  and 
(assuming  that  [ (0)=0)  find  the  first-order  results: 


A„  (G) 


a 


=  1 


1  C(3-G*)  [c^.A„(G‘)  -  agg.A_L(S*  )  ] 


-4 . 

G' 


CL-^ 

GG 

c-^ 

GG 


and 

Aj_(3) 


=  r  n3-3’)  [a33'A,|(3')  +  bgg.A^ca')] 


G' 


5 


(lla) 


(lib) 


Here  the  prime  on  the  summation  indicates  that  GVG,  and  we  have 
employed  the  notation  A||  (G)  for  A||(w,i^g)  (similarly  for  Aj^)  and 
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have  used  the  following  definitions: 


/3^  =  a(w,Kg)  +  a^(w,k||) 

(12a) 

=  a(w,K^)  -  a^(o,K^.) 

(12b) 

^GG'  "  ^3*  (^G  ^  ^G'  ) 

(12c) 

/\  A 

•  ^3- 

(12d) 

X  V  .  .  1  G  G 

^GG'  “  a(w,K^)a^(cj,K^.  ) 

(12e) 

< 

< 

®G  ^G  ^11  a(w,K^)a^((o,k||) 

(12f) 

r  il/2 

a  (w  ,Kg)  =  £  (w) -2  -  Kg 

(12g) 

As  in  Ref. 9,  we  assume  that  the  amplitudes  of  the  two 
components  of  the  specular  beam,  A||  (0)  and  Aj^(O),  can  dominate  the 
diffracted  orders.  Unlike  the  case  considered  in  Ref. 9,  we  now 
consider  the  possibility  of  having  four,  not  two,  resonant 
evanescent  waves.  We  may  have  four  wavevectors  in  the  surface 
plane,  corresponding  to  four  different  reciprocal  lattice  vectors, 
which  satisfy  the  resonance  condition: 


r 


sp 


for  r=l , . . . , 4 


(13) 


where  ^  surface-polariton  wavevector  corresponding  to 

the  frequency  of  the  incident  light.  To  zeroth-order  in  the 

roughness  (the  flat  surface  limit),  the  locus  of  points  ^  (w) , 

sp 

for  a  given  frequency,  defines  a  circle.  Hence,  to  zeroth-order 

the  four  vectors  lie  on  a  circle,  at  the  four-fold  resonance. 

If  any  one,  or  all  four,  of  the  are  at  resonance,  then  for 

GIT 

that,  or  those,  wavevectors  we  have  to  zeroth-order: 


8 


0 


(14) 


a(<o,Kg  )  +  e  ((j)a^((j,Kg  )  = 
r  r 

We  are  therefore  going  to  assume  that  six  amplitudes  can  be 
larger  than  all  the  rest,  these  being  A||(0),  Aj^(O),  A||  (G^)  , 
A||(G2),  A||  (G^)  ,  and  A||(G^);  and  for  each  of  these  we  can  write  an 
explicit  equation,  using  Eg. (11),  with  the  six  important  terms 
separated  out  from  the  sum  on  the  right  hand  side.  For  A||(l)  this 
procedure  gives  us: 


a 


A„(l)  = 


i  [  I  f(l-p)c^pA||(p)  +  C(l)[c^oA,|(0)  -  a^j,A^(0)] 

p=2,3,4 

+  C(l)[e^B„  -  +  {  ^’  C  d’j  )  [c^jA„  ( j  ) -a^^A^  ( j  )  ] 

j 

-  y  C  (1-p)  a  A  (p) 

p=2,3,4 


To  simplify  the  notation  we  replace  the  variable  or  subscript  G^ 
everywhere  by  j  alone  (where  j=l,..,4  are  the  resonant  beams  and 
j=0  the  specular  beam) .  The  prime  on  the  summation  indicates  to 
sum  on  all  j?£0,l,2,3,4  .  The  nonresonant  terms  have  now  been 

segregated  in  the  brackets  on  the  right  hand  side.  The  equation 
for  A||  (2)  is  obtained  from  Eq.  (15)  by  simply  interchanging  the  1 
and  2;  and  similarly  for  A||  (3)  and  A||(4). 

The  specular  amplitudes  are: 


“  “oo  1  ■ 

r=l 


a 


°°  R 


+  la 


00 


{Y  C(-j)[a(,jA„(j)+b„.A^ 


(j)]  +  }  C(-r)h^^Aj_(r)j  (16) 


r=l 
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(0,  =  C(-r)c„^A„(r, 


r=l 


+  !• 


i—  ly 

‘^oo  M 


f(-j)  [CojA||(j)-a^jAi(j)]  -}  C(-r)aj,^A^(r)}  . 


(17) 


r=l 


The  amplitucie  of  a  given  nonresonant  wave  can  also  be 

expressed  with  Eq. (11) ,  but  now  we  drop  all  the  other  nonresonant 

terms  from  the  sum  on  the  right  hand  side: 

4 

A_L(j)  =  i“jj{  J  C  (j-P)ajpA||  (P)  +  f  (j)  [a^^Aii  (0)+bjQAj_(0)  ] 
p=l 


+  C(j)[a.QB„+b.^BJ  } 


(18) 


A.iCj)  = 


la  .  .  /  ^ 

u 

p=i 


C  (j-P)CjpA||  (p)  +  f  (j)  [c^qAii  (0)-ajQAj_(0)  ] 

+  C(j)  }  . 


(19) 


We  now  follow  the  development  given  in  Ref. (9) .  We 
substitute  Eqs.(l8)  and  (19)  for  the  nonresonant  amplitudes  back 
into  Eq.(15)  for  the  the  resonant  amplitude  A(l);  and  for  the 


denominator  c^^  we  use 


K, 


‘^ll  ^  a(l)a  (1) 


a(l)a^(l)+K^ 

a(l)ao(l) 


(20) 


The 


where  1^.  now  stands  for  and  a(j)  for  a((o,Kri.) 

J  <^1 

denominator  c^^  vanishes  to  zeroth-order  at  the  i^^-resonance . 
We  multiply  both  sides  of  the  resulting  equation  for  A..  (1)  by 


( 1) (1) -K^,  in  the  form 

a2(l)a^(l)-K^  =  [a(l)a^(l)-Kj]  [a  ( 1 )  ( 1) +kJ  ] 


(21) 


on  the  right  hand  side  and  as 
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2  o)^ 

a2(l)a2(l)-Kj  =  e  ^  •  (22) 

c  c 

on  the  left  hand  side.  Here  we  have 

‘Jj.  =  (— )c  ,  (23) 

which  is  the  frequency  of  a  surf ace-polariton  at  wavevector  ^ 


on  the  flat  surface.  Finally,  we  note  that  in  the  equation  which 

we  now  have  for  A||(l),  we  are  left  with  an  '^||(1)  '^^e  right 

"2 

hand  side  in  a  term  proportional  to  C  ,  which  arose  from  the 

nonresonant  terms.  We  bring  this  term  over  to  the  left  side,  and 

use  it  to  renormalize  the  flat-surface  polariton  frequency: 

4 


-  2  2 
w  =  w 
r  r 


2 

C  a  nr 
rr  r 


(l+e)(j^K  ^ 

\  /  j- 


C(r-p)C (p-rjOppa 


2 

rp 


p=l 

(PJ^r) 


+ 


Y  nr-j)C(j-r)a.  . 

j 


^jj 


for  r=l,2,3,  or  4. 

Here  we  employ  the  definition: 

=  a  (r)Q:^(r)  [a  (r)a^(r)  -  ] 

Next  we  go  back  to  Eqs. (18)  and  (19)  for  the  nonresonant 
terms  and  substitute  them  into  Eqs.  (16)  and  (17)  for  Aj^(O) 


(24) 


(25) 


and 


A„ (0)  . 

At  this  point  we  have  six  linear  algebraic  equations  for 
A||  (0)  ,  Aj^(O),  and  A||  (r)  [r=l,...,4]  in  terms  of  B||  and  Bj^  ; 

expressed  in  matrix  form  these  equations  are: 
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S2 

S3 

S4 

Ml 

s 

A„(l) 

Si 

^23 

S4 

^2 

^2 

A||  (2) 

Si 

S2 

(n^-n^^) 

S4 

M3 

^3 

A„(3) 

Si 

S2 

S3 

M4 

S 

• 

A||(4) 

s 

s 

S 

S 

(i+p) 

Q 

A||  (0) 

^2 

^3 

^4 

s 

(l+T) 

Aj^(O) 

(U,-M^)B„  - 

NlS 

(s-s)s  - 

N2S 

(S-S)S  - 

N3S 

(s-s)s  - 

N4B1 

(W+V-P)B,|  - 

QBj^ 

a 

(X-S)B|,-  (T+- 

In  this  equation  Q  =  u/(2tTC/a)  =  a/X ;  =  w^/(2Trc/a);  while  the 

other  symbols,  namely  the  capital  letters  L  through  X,  are  defined 
in  the  appendix. 

The  problem  now  is  to  solve,  numerically,  this  set  of  six 
equations  for  the  six  unknown  coefficients,  the  A's.  Then  these 
six  coefficients  can  be  substituted  back  into  Eqs. (18)  and  (19)  to 
find  the  nonresonant  coefficients.  With  all  the  coefficients  thus 
at  hand,  we  can  then  determine  the  scattered  field  vector- 
amplitude,  ^^(StJ),  from  Eq.  (7)  —  for  each  diffracted  order  and 
for  each  evanescent  wave.  We  can  thus  obtain  the  reflectance,  or 
efficiency,  associated  with  each  of  the  grating's  diffracted 


12 


orders,  as  well  as  the  total  field  and  field-enhancement  near  the 
grating  surface. 

III.  Conclusion 

The  analytic  formulation  of  a  perturbation  theory  has  been 
presented  here,  in  complete  analogy  to  that  developed  in  Ref. 9  , 
in  order  to  treat  the  special  case,  of  normal  incidence  on  the 
bigrating  surface,  which  was  not  treated  in  Ref. 9  .  Numerical 
implementation  of  the  theory,  in  direct  analogy  to  that  in  Ref. 9 
is  now  being  carried  out.  This  implementation  will  include  a 
comparison  of  the  perturbation  theory  to  the  exact  ,i.e., 
non-perturbative,  theory,  in  order  to  test  the  limits  of  validity 
of  the  former.  When  the  conditions  at  normal  incidence  are  such 
as  to  excite  simultaneously  four  surface-polaritons ,  the 
absorptance  is  expected  to  be  enhanced  approximately  four-fold 
over  that  found  when  a  single  surface-polariton  is  resonantly 
excited. 

Again,  in  analogy  to  the  discussion  in  Ref. 9,  the  results 
presented  here  can  be  used  to  calculate  the  dispersion  relation 
for  surface-polaritons  on  bigratings.  We  take  the  matrix  on  the 
left  hand  side  of  Eq. (26) ,  cross  out  the  fifth  and  sixth  rows  and 
columns  (corresponding  to  the  specular  beams) ,  take  the 
determinant  of  the  resulting  4x4  matrix,  and  set  it  equal  to  zero. 
The  primed  summation  symbols  that  appear  in  the  definitions  of  the 
matrix  elements  must  first  be  redefined  to  reinclude  the  specular 
beams,  j=0.  Solution  of  the  determinantal  equation  by  fixing  iC|| 
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and  solving  for  complex  w  will  then  yield  both  the  dispersion 
and  the  life-time  of  the  surf ace-polariton .  The  present 
formulation,  unlike  that  in  Ref. 9,  will  enable  us  to  treat  those 
points  in  the  }C|| -plane  where  two  Brillouin  zone  boundary  lines 
intersect,  that  is,  where  there  are  four  coupled  surface 
polaritons.  Work  on  numerically  solving  the  dispersion  relation 
is  just  getting  started. 
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Appendix 

Eq. (26)  employs  the  following  definitions  for  the  matrix 
elements . 


L  =  h  if  (r-p)  c 
rp  r  [  '  '  rp 


-  Y  (f 


(A-l) 


j 


4 


pp  rp  pO 


p=l 

r) 


(A-2) 


j 
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N  =  -h 
r  r 


p=l 
(?^  r) 


“  y  f  (i^“j  )  C  ( j  )  ( f  •^•n+^  .b.^)a.. 

Z,  J/i.VJ/  \  jQ  j-j 


-a 


0  = 


00 


r  c 


00  '• 


ir(-r)Cor  +  J  C(p-r)C(-p)a  a 


p=l 

r) 


pp  pr  Op 


-  )  C  (j-r)f  (-j)  (f^.c.  -a  .a.^)Q;.. 
L  Oj  jr  rj  jO'  jj 

j 


00  ^ 


p=i  j 


00  I- 


np)c(-P)<.ppa^pbpc  ^  1'  c(j)C(-j)(foj-b^o)»jo- 


p=l 


=  '“OO 


iC(-r)aor  -  }  C (p-r) C (-p)a  a 


p=l 
{7^  r) 


PP  pr  pO 


-  y  C ( j “r) f (-j  )  (f  .  a^.+a.  b.„)a.. 
L  ^  '  jr  Oj  jr  3)0'' 

j 


=  “oo[  1 


C(P)f(-p)appap^bpjj  +  I  f(j)C(-j)  (fjo-bjo>®Oj“j 


p=l 

4 


*  “oo[  1 


C(P)C(-p)appbp^ 


p-1 


+  Y  nj)f(-j)(  ^  +b  ^  )a.. 
-1 


u 


=  h^[inr)g^  -  Y  C(r-j)C(j) 


T  .  C  . 

D  r: 


^  I  f<j)C(-j)T.C„. 


00  . 
3 


(A-3) 

(A-4) 

(A-5) 

j  j. 

(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A-10) 

(A-11) 
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w 


(A-12) 


X 


a 


c 


00^0 

oo^O 


“ooF 

j 


(A-13) 


In  the  above  equations  the  following  additional  definitions  have 

been  introduced: 
c  . 


f  . 

f . 
jm 


g  .  =  c  .  -  -  e  . 
^3  30  3 


h  =  - 
r 


nr  a  c 
r  rr 


6(0 


T  . 

3 


g  .a  .  ./c.  . 
3  33  33 


1 


1  2 


27rc/a, 


(A-14a) 

(A-14b) 

(A-15) 

(A-16) 

(A-17) 


The  terms  a...  b..,  c...  e.  were  defined  in  the  text.  Eg. (12). 
They  are  dimensionless,  as  are  f^j  and  g^  .  The  summations  2' 
are,  as  in  the  text,  over  all  j  except  j=0, 1, 2 , 3 , and  4  . 
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